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(57) Abstract 

A method, a map and an article of manufacture for the exploration of hydrocarbons. In one embodiment of the invention, the method 
comprises the steps oft accessing 3D seismic data; dividing the data into an array of relativey small three-dimensional cells; determining 
in each cell the semblance/similarity, the dip and dtp azimuth f the seismic traces c ntained therein; and displaying dip. dip azimuth 
and the semblance/similarity f each cell in the fonm of a two-dimensional map. In one embodiment, semblance/similarity is a function 
of time, the number of seismic traces within the cell, and the apparent dtp and apparent dip azimuth of the traces within the cell; the 
semblance/similarity of a cell is determined making a plurality of measurements f the semblance/similarity of the traces within the 
cell and selecting the largest of the measurements. In addition, the apparent dip and apparent dip azimuth, corresponding to the largest 
measurement of semblance/similarity in the cell, are deemed to be estimates f the true dip and one dip azimuth of the traces therein. A color 
map, characterized by hue, saturation and lightness, is used to depict semblance/similarity, true dip azimuth and true dip of each cell; true 
dip azimuth is mapped onto the hue scale, true dip is mapped onto the saturation scale, and the largest measurement of semblance/similarity 
is mapped onto the lightness scale of the color m^. 
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Method and Apparatus for Seismic Signal Processing and 

Exploration 

Technical Field 

This invention relates to the general subject of seismic exploration and. in 
particular, to methods and devices for identifying structural and stratigraphic features 
in three dimensions. 

Cross^Referencft 

This patent application is a continuation in part of a provisional patent 
application filed on October 6. 1995, and having a serial number of 005,032 and a 
U.S. Patent Application to Bahorich and Farmer, having a serial number of 353,934 
and a filing date of December 12.1994. 

Background Of The Invention 

In seismic exploration, seismic data is acquired along lines (see lines 10 and 
11 of FIG. 1) that consist of geophone arrays onshore or hydrophone streamer 
traverses offshore. Geophones and hydrophones act as sensors to receive energy 
that is transmitted into the ground and reflected back to the surface from subsurface 
rock interfaces. Energy is often provided onshore by Vibroseis® vehicles which 
transmit pulses by shaking the ground at pre-determined intervals and frequencies on 
the surface. Offshore, airgun sources are usually often used. Subtle changes in the 
energy returned to surface often reflect variations in the stratigraphic, structural and 
fluid contents of the reservoirs. 

In performing three-dimensional (3D) seismic exploration, the principle is 
similar; however, lines and arrays are more closely spaced to provide more detailed 
subsurface coverage. With this high density cov rage, xtrem ly larg volumes of 
digital data need to be record d. stored and processed before final interpretation can 
be made. Proc ssing requires extensive comput r resources and complex software 
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to enhanc the signal received from the subsurface and to mute accompanying noise 
which masks the signal. 

After the data is processed, geophysical personnel assemble and interpret th 
3D seismic information in the form of a 3D data cube (See FIG. 2) which effectively 
represents a display of subsurface features. Using this data cube, information can be 
displayed in various forms. Horizontal time slice maps can be made at selected 
depths (See FIG 3). Using a computer workstation, an interpreter can also slice 
through the field to investigate reservoir issues at different seismic horizons. Vertical 
slices or cross-sections can also be made in any direction using seismic or well data. 
Seismic picks of reflectors can be contoured, thereby generating a time horizon map. 
Time horizon maps can be converted to depth to provide a true scale structural 
interpretation at a specific level. 

Seismic data has been traditionally acquired and processed for the purpose of 
imaging seismic reflections for structural and stratigraphic interpretation. However, 
changes in stratigraphy are often difficutt to detect on traditional seismic displays du 
to the limited amount of information that stratigraphic features present in a cross- 
section view. While working with both time slices and cross-sections provides an 
opportunity to see a much larger portion of faults, it is difficult to identify fault surfaces 
within a 3D volume where no fault reflections have been recorded. 

Coherence is one measure of seismic trace similarity or dissimilarity. The 
more two seismic traces increase in coherence, the more they are alike. Assigning a 
coherence measure on a scale from zero to one. "0" indicates the greatest lack of 
similarity, while a value of "1" indicates total or complete similarity (i.e.. two identical, 
perhaps time-shifted, traces). Coherence for more than two traces may be defined in 
a similar way. 

One method for computing coherence was disclosed in U.S. Patent 
to Bahorich and Farmer (assigned to Amoco Corporation) 



having a serial number of 353,934 and a filing date of December 12, 1994. Unlike th 
shaded r li f m thods that allow 3D visualization of faults, channels, slumps, and 
other sedim ntary features from picked horizons, the coherency process devised by 
Bahorich and Farmer operates on the seismic data itself. \A/hen there is a sufficient 
Chang in acoustic impedance, the 3D seismic coherency cube dev loped bv 
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BahorichandFarm rcanbe xtremely ffecth/e in delin ating seismic faults. It is also 
quite ffectiv in highlighting sutrtle chang s in stratigraphy (e.g.. 3D images of 
meandering distributory channels. po.nt bars, canyons, slumps and tidal drainage 
patterns). 

Although the process invented by Bahorich and Famier has been very 
successful, it has some limitations. An inherent assumption of the Bahorich invention 
« the assumption of zero mean seismic signals. This is approximately true when the 
correlation window exceeds the length of a seismic wavelet. For seismic data 
containing a lOHz component of energy, this requires a rather long I00ms window 
whKrh can mix stratigraphy associated with both deeper and shallower time horizons 
Shortening the window (e.g.. to 32 ms) results in higher vertical resolution, but often at 
the expense of increased artifacts due to the seismic wavelet. Unfortunately, a more 
ngorous. non-zero mean running window cross correlation process is an order of 
magnmide more computationally expensive. Moreover, if seismic data is 
contaminated by coherent noise, estimates of apparent dip using only two traces will 
be relatively noisy. 



Thus, there is a need for methods and apparatus that would overcome the 
short-comings of the prior art. In particular, improved resolution and computational 
speed are desirable. In addition, it would be highly desirable to improve estimates of 
20 dip in the presence of coherent noise. 

Summary rtfj^i^ tnYmtinn 
In accordance with the present invention, a method and an article of 
manufacture is disclosed for locating subterranean features, faults, and contours In 
one embodiment of the invention, the method comprises the steps of: accessing 3D 
seismic data covering a presJetemiined volume of the ea^th. dividing the volume into 
an array of relatively small three-dimensional cells, wherein each of said cells is 
Characterized by at least five laterally separated and generally vertical seismic traces 
located therein: detemiining in each cell the semblance/similarity of the traces relative 
to two pre-detem^ined directions, and displaying the semblance/similarity of each cell 
in the fom, a two-dim nsional map. In on mbodiment. s mblance/similarity is a 
function of time, the number of seismic traces within the c II. and the apparent dip and 
apparent dip azimuth of the traces within the cell, the semblance/similarity of a cell is 
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determined by making a plurality of measurements of the s mblance/sirhilarity of the 
traces within the cell and sel cting the largest of the m asurements. In addition, the 
apparent dip and apparent dip azimuth, corresponding to the largest measurement of 
semblance/similarity in the cell, are deemed to be estimates of the true dip and true 
5 dip azimuth of the traces therein. Finally, a color map. characterized by hue. 
saturation and lightness, is used to depict semblance/similarity, true dip azimuth and 
true dip of each cell; in particular, true dip azimuth is mapped onto the hue scale, true 
dip is mapped onto the saturation scale, and the largest measurement of 
semblance/similarity is mapped onto the lightness scale of the color map. 

10 In another embodiment of the invention, an article of manufacture is disclosed 

that comprises a medium that is readable by a computer and that canies instructions 
for the computer to perform a seismic exploration process. In one embodiment, the 
computer accesses 3D seismic data covering a pre-determined volume of the earth 
and the medium instructs the computer to: divide the volume into an an^y of relatively 
15 small three-dimensional ceils, wherein each cell is characterized by at least five 
laterally separated and generally vertical seismic traces located therein; determine in 
each cell the semblance/similarity of the traces relative to two pre-determined 
directions; and store the semblance/similarity of each cell for display in the form a two- 
dimensional map. In one embodiment, the instructions on the medium define 
semblance/similarity as a function of time, the number of seismic traces wrthin the cell, 
and the apparent dip and apparent dip azimuth of the traces within the cell; the 
semblance/similarity of a cell is determined by making a plurality of measurements of 
the semblance/similarity of the traces within the cell and by selecting the largest of the 
measurements. In addition, the apparent dip and apparent dip azimuth, 
corresponding to the largest measurement of semblance/similarity in the cell, are 
deemed to be estimates of the true dip and tnje dip azimuth of the traces therein. The 
computer comprises means for producing a color display that is characterized by hue. 
saturation and lightness; and the medium has instructions to map true dip azimuth 
onto a hue scale, true dip onto a saturation scale, and the largest measurement of 
30 semblance/similarity onto a lightness scale. 

The process of the invention is particulariy well suited for interpreting fault 
planes within a 3D seismic volume and for det cting subtle stratigraphic features in 
3D. This is b cause seismic traces cut by a fault line generally have a different 
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seismic character than traces on eith r side of the fault. Measuring multi-channel 
coh rence or trace similarity along a time slice r v als lin aments of low coherenc 
along these fault lines. Such measures can reveal critical subsurface details that are 
not readily apparent on traditional seismic sections. Also by calculating trace similarity 
along a series of time slices, these fault lineaments identify fault planes or surfaces. 

The process of the invention presents a multitrace semblance method that is 
generally more robust in noisy environments than a three trace cross conflation method 
for estimating seismic coherency. In addition, the semblance process presented in this 
patent application provides: 

higher vertical resolution for good quality data than that of a three trace 
cross correlation measurement of seismic coherency: 

the ability to map the 3D solid angle (dip/azimuth) of coherent events; 

the ability to generalize the concept of complex "trace" attributes to one of 
complex "reflector attributes: and 

by combining these enhanced complex trace attributes with coherency 
and solid angle, the basis of quantitative 3D seismic stratigraphy data attributes 
that are amenable to geostatistical analysis methods. 

Moreover, seismic coherency versus dip maps of picked horizons allow 
analysis of: 

the structural and stratigraphic framework before detailed picking starts: 

structural and stratigraphic features of the entire data volume, including 
zones that are shallower, deeper, and adjacent to the primary zone of 
interest: 

subtle features that are not respresentable by picks on peaks and 
troughs: and 

features internal to the top and bottom of formation or s quence 
boundary picks. 
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Coupled With coherency, data cubes of th solid angi dip of coherent seismic 
reflection v nts allow one to quickly see structural as w II as stratigraphic 
relationships (such as onlap and offlap) between the seismic data and interpreted 
sequence boundaries. 

Numerous other advantages and features of the present invention will become 
readily apparent from the following detailed description of the invention, the 
embodiments described therein, from the claims, and from the accompanying 
drawings. 

Brief Desrriptjf^n pf the Dratyinqs 

FIG. 1 is a schematic diagram showing an an^ngement of geophones to 
obtain 3D seismic data from the earth s subsurface for processing in accordance with 
the present invention: 

FIG. 2 is a pictorial representation of the infonnation obtained from the data 
acquired using the arrangement of FIG. 1; 

FIG. 3 is a pictorial representation of a horizontal time slice (t=1200ms) of 3p 
seismic data processed in accordance with the prior art; 

FIG'S. 4A through 4H illustrate various analysis windows (computational stars) 
that may be used in running window analysis of seismic coherence, dip and dip 
azimuth: 

FIG. 5 is a pictorial representation of the process of the invention using an 
elliptical window centered about an analysis point; 



FiG s. 6A and 6B are examples of a rectangular dip/azimuth tessellation useful 
when analyzing a survey having strikes and dips parallel to the acquisition axes, and 
when illuminating faults cutting perpendicular to a dominant reflector strike and dip (po. 
25 qo); 

FIG'S. 7A through 7C ar pictonal repr s ntations of thr e t sselations of soUd 
angle dip/azimuth space; 
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FIG's, 8A through 8D depict the mapping of 3D seismic attributes (<J).c.d) to 3D 
color space (H.L.S); 

FIG. 9 shows four surfaces through the color hemisphere of FIG. 8A for four 
values of coherence; 

5 FIG'S. 10A through 10C depict ordinary vertical slices of the seismic data of 

FIG. 3; 

FIG'S. 11A through 11C depict the seismic attributes, dip, dip azimuth and 
coherency obtained by applying the process of the invention, to data corresponding to 
that of FIG.'s IDA through 10C; 

10 FIG'S. 12A and 12B are time slices (t=1200ms and t=1600ms) through the dip 

azimuth cube giving rise to FIG's. 11A and 11B; 

FIG'S. 13A and 13B are gray scale displays of coherency; 

FIG'S. 14A through 14C depict coherency slices corresponding to the data of 
FIG'S. 10A through IOC; 

15 FIG'S. 15A and 15B depict the results of applying a semblance algorithm and 

applying dip/azimuth algorithm in accordance with the present invention; and 

FIG'S. 16A and 16B are schematic diagrams depicting the processing flow of 
the steps performed in one embodiment of the invention. 

Detailed Description 

20 While this invention is susceptible of embodiment in many different forms, 

there is shown in the drawings, and will herein be described in detail, specific 
embodiments of the invention. It should be understood, however, that the present 
disclosure is to be considered an exemplification of the principles of the invention and 
is not intended to limit the invention to any specific embodiment or algorithm descnbed 

25 her in. 

Th first step of the process (See FIG. 16A) is to obtain a set of seismic data in 
the form of seismic signal traces distributed over a thr e dimensional volume of the 
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earth. Methods by which such data is obtained and reduced to digital form for 
processing as 3D seismic data are Icnown to those skilled in th art. 

The Semblance Process 

The next step is to generate a "coherence cube." This is done by applying a 
5 multi-trace semblance algorithm to the 3D seismic data. This algorithm may take 
many fomis. Whatever its form, its function is to compare the similarity of nearby 
regions of seismic data within the 3D seismic volume. This value (or attribute) serves 
as a rather robust estimate of signal discontinuity within geologic fomiations. as well 
as signal discontinuities across faults and erosional unconformities. 

10 We define an analysis grid (or computational star) to be either an elliptical r 

rectangular pattem of "J" traces centered about a given output trace (See FIG s. 4A 
through 4H). 

In the drawings, "X" denotes the center of the analysis window while "O" denotes 
additional traces used in the semblance calculation. Minimum size circular and 

15 rectangular windows used to analyze data with equal trace spacings (ax = Ay) are 
shown in FIG's. 4A and 4D. Minimum circular and rectangular windows used to analyze 
data ¥wth trace spacing in the cross-line/strike (y) direction twice that in the in-line/dip (x) 
direction (Ay « 2Ax) are shown in FIG s. 4B and 4E. Such nonequal spacings are 
commonly used to exploit the stower change of geology in the strike direction. Larger 

20 analysis windows used for greater resolution of reflector dip and azimuth, or to increase 
signal to noise ratio in poor data areas, are shown in FIG s. AC and 4F. 

Elliptical and rectangular analysis windows centered about an analysis point 
defined by a major axis, a minor axis, and the azimuth of major axis are shown in FIG s. 
4G and 4H. The acquisition (x.y) axes are retated by ^ degrees from the North-East (x'. 
25 y ) axes. Such assymmetric windows are useful in fracture detection. 

If we center the (x. y) axis about the center of an analysis window containing J 
seismic traces. Uj(t. Xj. yp. we define the semblance <t(t. p. q) to be: 
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c{T, p, q) = 



- (px, + qyj), X3, yj 

U°' 



- (p^i qyjX x^, y, ]}' 



(1) 



where the triple (t. p. q) defines a local planar event at time x. and p and q are the 
apparent dips in the x and y directions measured In ms/m. Since, p = d sin 4( and q = d 
cos ^. where d is the true dip and ^ is the dip azimuth, it follows that: 

5 Uj (t. p. q, X. y) * Uj [x - d(x sin «^ + y cos x. yj. 

Those skilled in the art will recognize that, in the denominator of equation (1), J 
serves as a normalization factor. The numerator represents the average energy and the 
summation term in the denominator represents the total energy of the traces. In effect, 
equation (1) is representative of a ratio of coherent and incoherent energy. 

10 The otqective is to perform a simultaneous 2D search (See FIG. 5) over apparent 

dips (p, q) in the in-line and cross-line directions. However, the semblance estimate 
given by equation (1) will be unstable for small but coherent values of seismic events, 
such as might occur if we were to sum along the zero crossings of a plane coherent 
wavelet. To avoid this, we estimate the coherence c (x, p, q) at time x and apparent dips 

15 (p, q) to be the average semblance over a time window (or vertical analysis window of 
height 2w ms of half length K = w/At samples): 



IC--K 



XuCt + kAl- (px, ♦ qyj),x^. yJ 



c(t.p.q) — ■ (2) 

JZZMt + kAt - (px, . qy.), Xj. y.]}^ 

In general, we do not know but wish to estimate that value of (p, q) associated 
with th local dip and azimuth of a hypoth tical 30 refl ction ev nt. 

20 In one embodiment of the process of the invention, we stimate (p. q) through a 

brute force search over all possible appar nt dips (S e FIG s. 6A and 6B). We assume 
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that the interpreter is able to estimate th maximum true dip. d^^, (measured in ms/m) 
from conventional s ismic displays of the data (e.g.. vertical data slices), thereby limiting 
the dips to be: 



*^ ^ "max- 



5 If 3W and )ffnax are the half width and half length of a rectangular analysis 

window, and if is the highest temporal frequency component contained in the 
seismic data, then the Nyquist criterion of sampling the data at two points per perioct 
restricts the apparent dip increments. Ap and Aq. to: 

x««^p ^ 1/(2f„«,). and y,^ Aq i V(2f^). 

10 It should be noted that the Nyquist criterion is valid for linear operations on the seismic 
data: and that equation (2) is nonlinear. In practice, we have found it necessary to limit 
Ap and Aq to half that required by the Nyquist sampling criterion to obtain an accurat 
semblance for a coherent dipping event. 

Thus, our search for an estimate of the apparent dip (p, q) of a seismic reflector 

15 is reduced to the calculation of semblance c(p,. q,„) over np * n, discrete apparent dip 
pairs 

(Pi. qm) where; 

np = (2d^Ap)+l.and 
nq = (2d^x/Aq)+1. 

20 The apparent dip pair (p,. q^) is deemed to be an estimate of the reflector 

apparent dips when: 

c(p.q) s c(p,. q„) (3) 

for all -Hp s 1 ^ ♦np. -n, s m s +nq. 
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The estimated apparent dips ( p . q ) are relat d to the estimated true dip d and 
dip azimuth ^ by the simple geometric relationships: 

p = d sin ^ ; and 

q s d cos ^ , 

5 where d is measured in ms/m and the angle ^ is measured clockwise from the positive 
x' (or North) axis. A simple coordinate rotation by angle is necessary when the in-line 
acquisition direction x is not aligned with the N*S (x') axis (See FIG. 4G). 

SouD Angle Discretization and Display 

Optimal angular discretization is important for two reasons: minimization of 
10 computational cost, and limitation on the number of colors that can be displayed using 
commercial interpretation workstation software (e.g., cunrently 64 with Landmark's 
"Seisworfcs" and 32 with Geoquest's "lESX" systems). 

FIG. 7A shows the discretization of apparent dip using equal increments Ap and 
Aq in a rectangular grid of 69 angles, FIG. 7B shows the discretization using equal 

15 increments Ad and A^ in a radial grid of 97 angles. Cleariy. we do not wish to sampi 
the dip d = 0 ms/m for ten different azimuths. The "Chinese Checker" tessellation of 
FIG. 7C more closely represents an equal and therefore more economic sampling of th 
(d. 4») surface with a minimum numt)er of points (i.e., 61 angles). Each tesseiation of 
FIG'S. 7A and 7C represents an approximately equal patch of solid angle AO. For the 

20 angular discretization shown in FIG. 7C and for a circular analysis radius. a« the 
incremental dip Ad is chosen to be: 

aAd^ :^ (4) 

xl max 

Display 

While it is possible to independently map s mblance. dip. and azimuth, it is clear 
25 that th latt r two attributes are coupled to each other. Furthermor , the confid nee we 
have in these estimates is proportional to the coher ncy/semblance. Oth rs (S e U.S. 
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Patent No. 4,970,699 to Bucher et al. and assign d to Amoco Corporation. "Method for 
Color Mapping G ophysical Data") have shown that the color HLS (hue. lightness, 
saturation) model can be quite effective in displaying multicomponent seismic attributes 
(Also see Foley, J.D. and Van Dam. A., 1981, Fundamentals of interactive Graphics 
5 Addison-Wesley, Reading. MA). 

Refering to FIG*s. 8A through 8D, in this scheme, we directly map azimuth. ^. 
onto the hue axis H: 

where both H (commonly known as the "color wheel") and <|> vary t>etween -ISO and 
10 +180 degrees (See FIG. 8B). Blue con-esponds to North, salmon to East, yellow to 
South, and forest green to West azimuth. Azimuths corresponding to zero dip are 
arbitrarily assigned a value of 0 degrees (North) and are thus plotted as blue. 

Next, we map (See FIG. 8C) average semblance/coherency c. onto the lightness 

axis L: 

15 L = ac. 

where 

O^L^IOO. 
O^c^ 1.0. and 

a is a scale constant less than 100. since changes in hue and saturation near L=0 
20 (black) and L=100 (white) are difficult to distinguish. White, or L = 100. corresponds to 
high semblance or c = 1. while black, or L=0, corresponds to low semblance, c = 0. 
Intermediate semblances correspond to intermediate shades of gray, (such as silver, 
gray and charcoal gray). Lightness (sometimes referred to as "brightness") expresses 
the amount of illumination. It represents a gray scale ranging from black to white. 

25 Finally, we map dip d onto the saturation axis S: 
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The saturation (S = 0) and hue chosen are arbitrary; we could just as asily have 
displayed this attribute for a value of (H = 0, S= 100) giving us semblanc display d as 
white, pastel blue, pure blue, midnight blue and black. Saturation expresses the lack of 
dilution of a color by white light. A fully saturated color has no white added; adding white 
5 "washes out' the color without changing its hue. (See FIG. 8D). 

FIG. 9 illustrates four constant surfaces through the 3D (H.L.S) color hemisphere 
of c. d) shown in FIG. 8A, corresponding to c = 1.00, c = 0.75. c = 0.50 and c = 0.00. 

Appendix 1 describes the color scheme in greater detail. Some advantages of 
the HLS color model are: azimuth is cyclic and maps neatly to the cydic color wheel 
10 (hue); the azimuths corresponding to d • 0 are meaningless; all azimuths converge 
smoothly to gray for shallow dips; and lower confidence in estimating dip and azimuth in 
zones of weak, low semblance (such as across faults) is indicated by darker colors. 

Implementation of Mathematical Process 

Landmark and GeoQuest interpretive workstations (See FIG. 16B). for example. 
15 can t>e used to view and interpret faults and stratigraphic features by loading the 
processed data as a seismic volume. Visualization software (e.g., Landmark's SeisCube 
software) may be employed to rapidly slice through the seismic volume to aid in 
understanding complex fault relationships. 

Computer Program 

20 A FORTRAN 77 program was written to perform the calculations and provide the 

information for the displays previously descrit>ed. Additional details are given in 
Appendix 2. Each trace Umn is accessed by its in-line and cross-line indices. M and N. 
The user specifies a rectangular or an elliptical spatial analysis window or cell about 
each point/trace in the input data set (See FIG. 4G). The major and minor axis of this 

25 analysis window, a and b are given by a - aplength and b = apwidth. The orientation 
or azimuth of the major axis is given by = apazim. A rectangular analysis window 
(FIG. 4H) is indicated by specifying -R on the command line. The 2J indices relative to 
the center of this analysis window (and corresponding to the traces that fall within this 
window) are tabulated as a simpi list, with m(j) and n(i) indicating the trace index 

30 (relative to the analysis trace Umn) in the x and y directions, respectively. The program 
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performs a simultan ous 2D search over apparent dips (p.q) in the in-line and cross-line 
directions, where (p^ + q^)^'^ < -i-smax. The increments dp and dq are chosen such that 
the data are sampled at four points per period =1/(fref) at the edge of the analysis 
window. For interpretation, it may be convenient to express each apparent dip pair (p.q) 
5 in spherical coordinates as a true (time or depth) dip d and dip azimuth ^. 

The data in the analysis window are interpolated to the fractional time, x - px - qy. 
for each trial dip and azimuth (See FIG. 5), in essence, "flattening" the data. The 
semblance for this trial dip at the analysis point is defined to be the semblance of these 
flattened traces in the analysis window. 

10 For time domain data, we flatten the jth trace about the analysis point (M.N) by: 

Uf{x. P. X. y)=u( T - (px+qy) ] =u(t - d(x sin* -►y cos*)]. 

where x and y are distances measured from the center of the analysis window. This 
may be expressed: 

15 where ^ and Ay are the irvline and cross-line trace spacings. 

For depth domain data we flatten the jth trace using: 

"K^. p. q. X. y)=uK - (px+qy) 1 =u[ ^ - d(x sintt* +y cos*)]. 
The semblance is then calculated for all subsequent dips and azimuths using: 



■ > V 



(5) 

20 As in velocity analysis, the semblance for each dip, azimuth and analysis point 

are smoothed by forming a running window time integration over the partial sums from - 
K to +K wh re K = aph ight/dt We therefore define the coh rence. c(T.p,q) to be: 
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+Kf +J 

_I^I[.f(/,p,^,x,>.)] 




(6) 



10 



15 



20 



That dip and azimuth pair n « (d. ^) which has the maximum (running window 
integrated) coherency c is taken to be an estimate of the coherency, c . dip and azimuth 
( d . ^ ) for the analysis point. 



FIG'S. 11 A through 11C are displays of the 3D seismic attributes c. d) 
corresponding to FIG's. 10A through IOC using the semblance based coherency 
algorithm expressed by equation (6). and the color display technique depicted in FIG's. 8 
and 9. The input data were temporarily sampled at 4 ms, have an in-line trace spacing 
of Ax = 12.5 m, and have a cross-line trace spacing of Ay = 25 m, with the in-line 
acquisition oriented along a N-S axis. For FIG's. 11A through 11C, a circular analysis 
window or cell of a - b - 60 m was used (See FIG. 4A). so as to include a total of 1 1 
traces in the calculation. The maximum search dip (See FIG. 7C) was dm^i = 0.25 
ms/m, giving rise to 61 search angles. The temporal integration time used was w = 16 
ms, or K s 4, thereby averaging the semblance calculation over 9 samples. 

In FIG'S. 10A and 10B lines AA' and BB* were chosen as S to N and W to E 
vertical slices through the center of a salt dome. Line CC is an offset S to N line and 
illustrates the apF>earance of radial faults on a vertical slice. In FIG's. 1 1 A through 1 1C. 
the interior of the salt dome is represented by dark colors, corresponding to an area of 
generally low coherency. Low areas of coherency correspond to the radial faults seen 
on line CC*. Coherent, flat dips are represented as light gray and dominate the section 
away from the salt dome, in particular line CC. The blue color on the north side of the 
salt dome (seen on N-S line AA ) corresponds to sediments dipping steeply (d = dmax) to 
the North. Thes dips l>ecom progressively shallower away from the salt dome, and 
are thus displayed first as blu (saturation. S = 100.0), cadet blue (S * 0.75) and steel 
blue (S = 0.50). before they flatten and are displayed as gray (S = 0.0). The yellow color 
on the south side of the salt dome (s en on line AA*) corr spends to s diments dipping 
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Steeply to the South. Th salmon color on the East flank of the salt dome (shown on the 
E-W line BB') con- sponds to s diments dipping steeply to the East. These dips also 
become progressively shallow away from the salt dome, and are displayed first as 
salmon (S = 100.00). through sienna (S = 50.0), and finally to gray, corresponding to flat 
5 dip. Finally, the forest green color on the West flank of the salt dome (shown on line 
AA') corresponds to sediments dipping steeply to the West. These dips also flatten 
away from the salt dome and are displayed using the colors shown on the West part of 
the legend shown in FIG. 9. N-S line CC is not aligned radially with the salt dome. 
Thus, out-of-the-plane rotation of different fault blocks are depicted, with the green block 
10 corresponding to dips to the SW and the cyan block with dips to the NW. 

Since these 30 attributes were calculated for every point on the input seismic 
volume, they can be displayed as horizontal attribute time slices (See F4G's. 12A and 
12B); these correspond to a time slice of the unprocessed seismic data. The interior of 
the salt dome, as well as the radial faults are displayed as dark colors, corresponding to 

15 incoherent zones of the data. Because of the neariy radial symmetry of the salt diapir at 
t = 1,200 ms (See FIG. 12A). the dipping sediments that flank the diapir also radiate 
outward in an azimuthally simple fashion such that their azimuths correspond quite 
closely to the color legend on the left side of FIG. 9. This pattem is somewhat less 
symmetric at t = 1.600 ms (See FIG. 12B). where there are shallower dips to the South 

20 than to the North. In addition, internal blocks of coherent data can be seen within th 
salt dome. 

The color legend displayed in FIG. 9 allows for only four "buckets" of coherency, 
in order to examine the coherency in greater detail, it can be plotted as a single attribute. 
This is shown in FIG'S. 13A and 138 where all 184 colors are applied to the simple gray 

25 scale shown of FIG. 8C. In this display, maximum coherency (c = 1.0) is rendered as 
white: minimum coherency (c = 0.0) is rendered as black. While the interior of the salt 
diapir is shown as a highly incoherent zone, this display t>etter shows subtle details in 
the radial faults patterns. In particular, faults emanating from the salt dome are shown, 
with some bifurcating as we move away. In addition to more continuous binning of the 

30 coherency attribute, part of this difference in perception is due to the fact that th human 
retina sees colors and black and whit using different (cone vs. rod) receptors. Ther is 
also a physiological difference in the ability to differentiate t>etween gre ns and btu s 
between male and f male populations. For this reason, male interpreters often prefer 
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the simple single attribute coh rency display shown in FIG's. 13A, 13B and FIG. ISA 
over the multiattribute (♦x.d) display shown in FIG's. 11A through 12B and FIG. 15B. In 
actuality, these displays are quite complimentary: the 3D component display being 
useful in recognizing the appearance of conflicting dips azimuths between adjacent 
5 rotated fault blocks; and the single component display being used to enhance the edge, 
or incoherent fault discontinuity, separating them. 
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Process Considerations 

Careful study of FIG s. 13A and 13B reveals a ring-like pattern of incoherent 
energy circumscribing the salt dome. To investigate the cause of these artifacts, vertical 
slices were taken through the single component coherency cube con-esponding to the 
5 seismic data in FIG s. 10A through 10C. This is shown in FIG s. 14A through 14C. The 
interior of the saK dome is cleariy incoherent. An incoherent submarine canyon feature 
(described by Nissen et al., "3D Seismic Coherency Techniques Applied to th 
Identification and Delineation of Slump Features". 1995 SEG Expanded Abstracts, 
pages 1532-1534) is shown to the north of the salt dome. If the seismic data shown in 

10 FIG'S. 10A through IOC were overiayed on the coherency sections shown in FIG s. 14A 
through 14C. one would see a close correspondence between areas of low coherency of 
FIG'S. 14A through 14C with zero crossings of the seismic reflection events in FIG's. 
IDA through 10C. This is easily understood if it is assumed that there is a fixed, but 
incoherent, level of seismic noise throughout the data. For analysis points where th 

15 apparent dtps are aligned with the peaks or troughs of strong amplitude seismic 
reflectors (such that the estimate of signal energy is high with respect to the incoherent 
noise), one can expect the signal-to-noise ratio to be high, giving rise to an estimate of 
high coherency. However, if. our analysis point is such that there are apparent dips 
aligned with the zero crossings of these same seismic reflectors, such that the signal is 

20 low with respect to our incoherent noise, one can expect the stgnal-to-noise ratio to be 
low, giving rise to a low estimate of coherency. 

We have found three methods for increasing the signal-to-noise ratio: the first 
more appropriate for structural analysis; the second more appropriate for stratigraphic 
analysis, and the third appropriate for both. 

25 For the case of steeply dipping (less than 45 degrees from the vertical) faults, the 

signal-to-noise ratio can be increased by simply increasing the size of our vertical 
analysis window w given in equation (2). Two effects will be observed. First, the 
structural leakage corresponding to the zero crossing points of the reflectors diminishes 
as vertical integration window size increases. S cond. since f w of the faults are truly 

30 vertical, the lateral r solution of the faults appears to decrease as the vertical window 
size increas s An analysis window of w = 16ms (which would ncompass a full cycle of 
the peak 30 Hz energy in the data) appears to be in good compromise. 
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The second method (equally appropriate for stratigraphic and structural analysis) 
of incr asing the signal-to-noise ratio, is to extract coherency along an interpreted 
stratigraphic horizon. If this stratigraphic horizon is associated with an extremum of the 
seismic data, such as a peak or trough, those data having only a relatively high signal- 

5 to-noise ratio are selectively displayed. Cleariy. extracting coherency data 
corresponding to a zero crossing would greatly exacerbate the coherency display. A 
more economic version of this approach is to first flatten the data along the horizon of 
interest and then calculate the seismic attributes only along the picked horizon. This 
approach is somewhat more sensitive to busts in automatic (and human!) pickers, since 

10 cycle skip glitches in the picking are somewhat random and therefore will almost always 
appear as incoherent. 

Shallow features (e.g.. shallow channels; shallow tidal channel features 
corresponding to reworked deltaic sands; and small en echelon faulting) do not exist for 
any distance above or below an interpreted stratigraphic horizon; therefore, the inclusion 

IS of any data from above or below the horizon in which they are located adds uncorrelated 
amplitude variations, thereby making these discontinuities look more coherent, and 
hence washed out. If the time samples above or below the interpreted horizon contain 
independent, perhaps strong amplitude discontinuities, these discontinuities will bleed 
into the analysis for large windows, giving a stratigraphic horizon containing featur s 

20 mixed from stratigraphic different horizons generated at different geologic times. 

The third method is a generalization of the original collection of seismic traces Uj 
to that of an analytic trace v, defined as. 

V,(t) s Vj(t) * iUi"(t) 

where u^"(t) is the quadrature, or Hilbert transform of Uj(t). and i denotes ^PT. The 
25 calculation of o(x.p,q) and c{t,p,q) is entirely analogous to equations (1) and (2), where 
we note that the definition of Vj^ is given by 

v,^sv,v,-s(Ui+iV|)(uj-iVj). 

The third method avoids numerical instabilities in the s mblance estimate of equation (1) 
at the "zero-crossings" of an otherwise strong reflector. 
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The Effect of the Horizontal Analysis Window 

By examining equation (2). it is clear that the computational cost of analysis 
increases linearly with the number of traces included in the analysis. However, by 
comparing a semblance based 11 -trace coherency time slice with those of a 3-trace 
cross correlation coherency time slice, (where each has an identical vertical analysis 
window of w = 32 ms) one is led to believe that adding more traces to the computation 
can increase the signal-to-noise ratio. In generat the signal to noise ratio increases as 
we increase the size of the analysis window. However, the overall coherency decreases 
somewhat (one sees less white), since the approximation of a possibly curving reflector 
by a constant (p.q) planar event breaks down as we increase the window size. In 
general, the signal-to-noise ratio of dip/azimuth estimates increases with the number of 
traces in the calculation, until a point is reached whereby the locally planar reflector 
approximation no longer holds. 

Conclusions 

15 The 3D semblance technique presented in this patent application provides an 

excellent measurement of seismic coherency. By using an arbitrary size analysis 
window, we are able to balance the conflicting requirements of maximizing lateral 
resolution and signal-to-noise ratio that is not possible when using a fixed three trace 
cross correlation technique. Accurate measurements of coherency can be achieved by 

20 using a short temporal (vertical) integration window that is on the order of the shortest 
period in the data, whereas a zero mean cross correlation technique preferably is used 
with an integration window that is greater than the longest period in the data. Thus, the 
semblance process results in less vertical smearing of geology than a cross con-elation 
process, even for large spatial analysis windows (See FIG's 15A and 15B). Equally 

25 important to the coherence estimate, the semblance process provides a direct means of 
estimating the 3D solid angle (dip and azimuth) of each reflector event. These solid 
angle maps may or may not be related to conventional time structure maps defining 
formation boundaries. Like the basic coherency process of Bahorich and Fanner (e.g.. 
cross correlation), estimation of the instantaneous dip/azimuth cube can be achieved 

30 prior to any int rpretation of the data for use in a gross overvi w of the geologic setting. 
In this reconnaissance mode, the coherency and instantaneous dip/azimuth cubes allow 
the user to pick key dip and strike lines crossing important structural or sedimentologic 
features very arty in the interpr tation phase of a project. In an interpretation mode, 

RECTIFIED SHEET (RULE 91) 
ISA/EP 

BNSDOCID: <WO 9713ie6AlJ_> 



wo 97/13166 



•21 • 



PCT/US96/15705 



10 



15 



these dips and azimuths may be related to fomiation and/or sequence boundaries, such 
that on can map progradation and transgr ssion patterns of the int mal structure in 3D 
Finally, having estimated the instantaneous dip and azimuth at every point in the data 
cube, one can apply conventional seismic trace attributes to locally planar reflectors, 
thereby greatly increasing signal-to-noise ratios. 

From the foregoing description, it will be observed that numerous variations, 
alternatives and modifications will be apparent to those skilled in the art. Accordingly, 
this description is to be construed as illustrative only and is for the purpose of teaching 
those skilled in the art the manner of carrying out the invention. Other algorithms may 
t>e used to measure the similarity of nearby regions of seismic data or to generate the 
"discontinuity cut>e/* Moreover, equivalent computations may be substituted for those 
illustrated and described. For example, instead of a search over apparant dips p and q. 
one could search over dip and azimuth (d, ^ ). The inverse of the computed semblance 
may t>e used so as to obtain a display analogous to the negative of a photograph. Also 
certain features of the invention may be used independently of other features of the 
invention. For example, after the solid angle (dip and azimuth) has been estimated, a 
smoother and more robust muititrace estimate of the conventional complex trace 
attributes (Taner, M. T., Koehler. F., and Sheriff. R. E.; 1979; "Complex Seismic Trace 
Analysis;" Geophvsies. 44. 1041-1063) may be obtained. Instead of calculating these 
attributes on a single trace, one can calculate attributes of the angle stack of traces 
within the analysis window. That is, one can calculate: 




V.('c,p.^) = tan 



-*{f/''(T,/7.^)/f/(T,p,9)}. 





and 



[t/(T,p,^)J^^[C/''(T,/7,9)]' 
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where 



5 U(T. p. q) is jXuC T - (px, + qy^), , y J 



10 



(See the numerator of equation 1); 



U"(t. p. q) is the Hilbert transform, or quadrature component of U(t. p. q); 

SiCt. P. Q) is the envelope, or instantaneous amplitude; 

U'i(t. P. q) is the instantaneous phase: 

fi(t. p, q) is the instantaneous frequency: and 

bi(t. p. q) is the instantaneous bandwidth (See Cohen, L.: 1993; 'Instantaneous 
Anything;- Proc. IEEE Int. Conf Anniist Snaftrh Signal ProM>««ing 4, 105-109). 



In addition to these "instantaneous* attributes, other attributes are suggested to 
characterize the signal within a given lobe of the trace envelope to be that of the 
attribute at the peak of the envelope x,. These include (See Bodine. J. H.; 1994; 
IS -Wavefomi Analysis with Seismic Attributes;" presented at the 54th Ann. Intl. Mtg. SEG. 
Atlanta. GA. USA): 
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the wavelet envelope: 

Br(x, p, q) = a, (xe, p. q). 

the wavelet phase: 

Vr(x» p. q) = Vi(Te, p. q). 

S the wavelet frequency: 

fr(x. p. q) = fi(ie, p. q). 

the wavelet bandwidth: 

br(x, p. q) = bi(xe, p, q). 

the zero phase component: 
10 U^(T , p,q) = cos[m; , (t . p, q)]U(T . p,q) + sin[ (t , p,q)]U" (x . p,q) 

the ninety degree phase component: 

U*(T,p.q) = -sin[v,(T,p,q)lU(T,p.q) +cos[v,(T,p,q)]U"(T,p,q) 

as well as skewness. rise time, and response length. Since mixing occurs along the 
true dip direction, slowly varying amplitude, phase, frequency, and bandwidth 

15 components of the event will be preserved. Moreover, the computation of 
coherency/semblance/similarity allows one to perform "texture analysis" of similar 
seismic regions. Texture analysis combined with "cluster analysis" leads to 
segmentation analysis. Among other things, this allows one to make geologic 
correlations and extrapolate the geological character of the subsurface. In addition, 

20 determination of the coherency may be used to impose q priori constraints for both 
post-stack and pre-stack s ismic inversion. Thus, it will be appr ciated that various 
modifications, alternatives, variations, and chang s may be mad without departing 
from the spirit and scope of the invention as defined in the append d claims. It is. of 
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course. intended to cover by the appended claims all such modifications involved 
within the scope of the claims. 
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Appendix 1 

MUL-nATTRiBUTE HLS CAUBRATION 

The hues are pure, or 100% saturated colors, and correspond to the following 1994 



non-toxic 96 crayon "Crayola" standard: 




direction 






N 


0 


Dlue 


NNE 


30 


plum 


ENE 


60 


magenta 


E 


90 


salmon 


ESE 


120 


red 


SSE 


150 


orange-red 


S 


180 


yellow 


SSW 


210 


Itme-green 


WSW 


240 


green 


W 


270 


forest-green 


WNW 


300 


cyan 


NNW 


330 


cerulean 


N 


360 


blue 


Partial 50% saturation corresponds to "dirtier* or "muddier* colors: 


diraetion 


^ (hue) 


Crayola Color 


N 


0 


cadet blue 


NE 


45 


fuscta 


E 


90 


maroon 


SE 


135 


sepia 


S 


180 


gold 


SW 


225 


olive 


W 


270 


sea green 


NW 


315 


steel blue 


N 


360 


cadet blue 


0% saturation corresponds to no color pigment: 




direction 


A (huet 


Crayola Color 


N 


0 


gray 


E 


90 


gray 


S 


180 


gray 


w 


270 


gray 


N 


360 


gray 



Low values of lightness correspond to "dark* colors: intermediate values of lightness 
correspond to "deep' colors, and high values of lightness correspond to "pastel' colors. 
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Appendix 2 

SYNOPSIS 

\ sembJd [ -^file^in ] [ Oftlej>ut] [ -Msfile^hls \ [ -t»tartcscarc ] { -tend tend ] 
[-ildmdx ] [-cidmdy l [ - aplangthapiengch ] [- apwidthapwidch] ( - apheightapiieighcl 
[ -apazimapaziin ] [ -ilazimxazim ] [ -clazimyaziin ] [ -dzdz ] [ •smaxsmax ] [ •pmtnpmiii ] 
5 ( -pmaxpinajc ] [ -qmingmin ][ -qmaxginajc ]( -thresh ch^-esh ] [-Inlfref] 
[ -startlifiascarciine] [ -andlineeiidiine] { -exppower] [-minj { -Int] [-R] 

DESCRIPTION 

semb3d reads in 3D seismic post stack time or depth data and generates semblance, dip 
and azimuth outputs. 

10 COMMAMD UNg ARGUMgMTS 

aainbSd gets ail its parameters from command line arguments. These arguments specify 
the input, output, spatial analysis window, and dip discretization parameters. The following 
command line arguments have been used in one embodiment of the invention. 

44file_in 

15 Enter the Input data set name or file immediately after typing -N. This input file should 
include the complete path name if the file resides in a different directory. Example: - 
N/export/data2/sanJuan/t:iine_stack tells the program to look for file ' time^scack* in 
directory 7export/data2/8anJuan'. For this program, the data is stored as a rectangular 
grid of regutariy binned data. The number of traces (denoted by lineheader word 'NumTrc') 

20 defines the number of traces in the *x* direction. The number of records (seismic lin s 
denoted by lineheader word 'NumRec*) defines the number of traces in the y direct! n. 
Missing data padded in with dead traces flagged by a dead trace header flag. 



-Ofile_out 

25 Enter the output multi-attribute data set name or file immediately after typing -O. Attributes 
will be output back to back, line by line. Without scaling the semblance c will rang 
t>etween 0.0 and 1.0. The values of dip will range between 0 and smax and will always be 
positive (pointing down). Units are in msec/m (msec/ft) for time data, or m/m (ft/ft) for depth 
data. The azimuth ^ is perpendicular to strike and points in the direction of maximum 

30 positive dip (pointing down). The values of azimuth will range between 0 and 360 degrees. 
Property defined, an output azimuth of 0 degrees corresponds to North, while an output 
azimuth of 90 degrees corresponds to East. The values of OMEGA s (d. 4>) can be chosen 
such that (when converted to an 8 bit integer) the left most 6 bits correspond to a valid 
Seisworks color table. This color table corresponds to the HLS color model previously 

35 descrik>ed and is generated using a program that maps the angles scanned into an HLS 
(hue. lightn ss. saturation) color map of OMEGA = (d, 4>). 
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-hls frie^hls 

Enter -his followed by the his table file name to output an ascii flat file containing the hue, 
lightness and saturation of each sample contained in the output. This file is input to a 
5 program to generate a RGB (red, green, blue) color lookup table needed for a proper 
display on certain workstations. 

-tstaittstart 

Errter -tstart followed by the beginning of the analysis window in msec. 

10 

-tendtend 

Enter -tend followed by the end of the analysis window in msec. The output record will be 
(tend - tstart) msec long. 

15 -ildmdx 

After -ildm enter the in-line distance measure (trace separation) in m (ft), 
•cldmdy 

After •cidm enter the cross-line distance measure (line separation) in m (fl). 

20 

-dzdz 

After -dz enter the vertical depth sample increment in m (fl). A value of dz > 0 indicates the 
data are in depth. 

25 -aplengthapleng^th 

After -aplength enter the half aperture length (in meters or feet) along the azimuth of the 
elliptical analysis window to be used. Increasing the analysis window by increasing 
aplength. apwidth will result in: (1) increased angular resolution. (2) decreased spatial 
resolution, (3) increased computational cost; and (4) decreased overall coherency (since the 
30 plane wave approximation is less valid. 
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-apwidthapwi d th 

After -apwidth enter the half half aperture width (in meters or feet) perpendicular to the 
azimuth of the elliptical analysis window to be used. 



5 -apheightapheight 

After -apheight enter the half length in milliseconds (or meters or feet) of the running time 
(depth) integration window applied over the semblance. Example = ±2 samples. Increasing 
the temporal Integration window apheight will result in: (1) a smoothed, less noisy respons . 
(2) decreased vertical resolution, and (3) no change in computational cost. 

10 

-apazimapaziii) 

After -apaatim enter the azimuth of the elliptical analysis window (with 0 being North and 90 
being East). 



15 •smaxsmax 

After -smax enter the maximum dip to be tested in msec/m (msecm) for time data, or in 
m/m (ft/ft) for depth data. This is recommended wtien there is no preferential strike 
direction in the data. This value can be read directly from a section display of the data. 
smax will be on the order of .30 msec/m (10 msecAfl) for time data. Increasing the value of 
20 smax beyond any true dips results in significantly increased computational cost for an 
identical result. 



-pminpinin 

After -pmin enter the minimum inline (increasing trace number) dip to be tested in msec/m 
25 (msec/ft) for time data, or in m/m (ft/ft) for depth data. This is recommended when there is 
a predominant strike direction parallel or perpendicular to the data acquisition lines. This 
value can be read directly from a section display of the data. 



•pmaxpniax 

30 After -pmax enter the maximum in-line (increasing trace number) dip to be tested in 
msec/m (msec/ft) for time data, or in m/m (ft/ft) for depth data. This is recommended wh n 
there is a predominant strike direction parallel or perpendicular to the data acquisition lines. 
This value can be read directly from a section display of the data. 
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-R 

Enter this command line argument to define a rectangular (2*aplength by 2*apwidth) vs. 
elliptical analysis window oriented along the azimuth axis. 

5 

•qmingmiii 

After -qmin enter the minimum cross*line (Increasing line number) dtp to be tested In 
msec/m (msec/ft) for time data, or in m/m (ft/ft) for depth data. This is recommended when 
there is a predominant strike direction parallel or perpendicular to the data acquisition lines. 
10 This value can be read directly from a section display of the data. 



-qmaxgmajc 

After -qmax enter the maximum cross-line (increasing line number) dip to be tested in 
msec/m (msec/ft) for time data, or in m/m (ft/ft) for depth data. This is recommended when 
15 there is a predominant strike direction parallel or perpendicular to the data acquisition lin s. 
This value can be read directly from a section display of the data. 



-thresh triiresh 

After -thresh enter the threshhold or cutoff semblance value, below which dip and 
20 azimuth are considered to be valid measures; below this value shades of gray will be 
displayed. Some display software limits the number of colors available for display. 



-freffref 

After -ffreff enter the reference frequency in cycles/sec (Hz) for time data, or in cycles/km 
25 (cycles/kft) used in detemiining the number of dips to be searched (e.g., fref = 60Hz for 
time data. 30 cycles/km for depth data). 



-ilazimiiazim 

After *ilazim enter the in-line azimuth (0 degrees being North. 90 degrees being East) that 
30 is the azimuth of increasing trace numt)er. This value is used to calibrate a solid angle 
output file, if used. 
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-clazimci a zim 

After -clazim enter the cross-line azimuth (0 degrees being North, 90 degrees being East) 
that is the azimuth of increasing line numbers. This value is used to calibrate the solid 
angle output file, if used. 

S 

•exppower 

After «exp enter the exponent to be applied for non*iinear scaling of the semblance. In 
general, most semblance/coherency values will be between 0.8 and 1.0. Scaling with 
power = 2.0 would map these values between .64 and 1.0, scaling with power = 4.0 would 
10 map these values between .41 and 1.0. and so forth. This is useful for loading data to an 
interpretive workstation. 



-startlinescartline 

After -startline enter the first output line to be generated. 

15 

•endlineendline 

After -endline enter the last output line to be generated. 



•min 

20 After -min enter this command line argument to extract the dip. azimuth, and semblance 
corresponding to the minimum semblance of the angles searched. (As a default, th 
program searches for the maximum semblance or coherency). 

-int 

25 Enter this command line argument to scale output such that it can be represented by an 8 
bit integer ranging between -128 and «M27. Useful for loading data to an interpretive 
workstation. 
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CLAIMS 

We claim: 

1 . A method for the exploration of hydrocartwns, comprising the steps of; 

(a) ot>taining a representation of a set of seismic traces distributed 
over a pre-determlned three-dimensional volume of the earth, said volume of the earth 
havmg subten^nean features characterized by dip and dip azimuth that are defined 
relative to a pre-defined dip azimuth measurement axis; 

(b) dividing said three-dimensional volume into at least one horizontal 
time layer, and dividing said time layer into a plurality of three-dimensional analysis ceils, 
wherein each analysis cell has two pre-determined. mutually perpendicular lateral 
dimensions and has portions of at least five laterally separated seismic traces locat d 
therein: 

(c) calculating, within each of said analysis cells, a plurality of 
measures of the semblance of said traces located therein, wherein each measure of 
semblance is at least a function of time, the number of seismic traces within said analysis 
cell, and the apparent dip and apparent dip azimuth of said traces within said analysis 
cell; 

(d) identifying, within each analysis cell, the largest of said calculated 
measures of semblance and defining the corresponding apparent dip and apparent dip 
azimuth to be an estimate of the true dip and an estimate of the true dip azimuth of th 
seismic traces within said analysis cell: and 

(e) forming, from all of said analysis cells, a seismic attritHJte display 
from said largest calculated measures of semblance and said corresponding estimates of 
the true dip and the true dip azimuth of the seismic traces within said time layer. 
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2. The method of Claim 1 , where step (e) is perfomied by fomning a color 
map that is charact rized by hue« saturation and lightness. 

wherein one of said estimates of true dip azimuth, said estimates of true dip. and 
said largest calculated measures of semblance is mapped onto one of a lightness scale. 
5 hue scale, and a saturation scale; 

wherein another of said estimates of true dip azimuth, said estimates of true dip. 
and said largest calculated measures of semblance is mapped onto another of said 
lightness scale, said hue scale, and said saturation scale; and 

wherein the remaining one of said estimates of true dip azimuth, said estimates of 
10 true dip. and said largest calculated measures of semblance is mapped onto the 
remaining one of said lightness scale, said hue scale, and said saturation scale. 

3. The method of Claim 2. where step (e) is performed by mapping said 
estimates of true dip azimuth onto said hue scale. 

4. The method of Claim 2. where step (e) is performed by mapping said 
15 estimates of true dip onto said saturation scale. 

5. The method of Claim 2. where step (e) is performed by mapping said 
largest calculated measures of semblance onto a lightness scale. 

6. . The method of Claim 1. where in performing step (c) each measure of 
semblance is at least a function of the energy of said traces; and wherein said energy of 

20 said traces is a function of time, the number of seismic traces within said analysis cell, 
and the apparent dip and apparent dip azimuth of said traces within said analysis cell. 

7. The method of Claim 6. wherein each measure of semblance is at least a 
function of: 

I J 
J-l J-» 
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where each analysis cell contains portions of at least J (J25) seismic traces, wh re x and 
y are distances measured from the center of the analysis cell, wh re p and q are the 
apparent dips in the x and y directions respectively, and where Uf (t. p. q. x. y) is a seismic 
trace within the analysis cell; and wherein the true dip d and dip azimuth (p are related 
to p and q by p = d sincp and q ^ d cos<p. 

8. The method of Claim 7. wherein each measure of semblance is a function 

of: 



(Z Uf(^P.q»Xj,yj)) 



J-i 



Z "f(^P.q»X3,yj)' 

9. The method of Claim 7. wherein each measure of semblance for each dip. 
10 dip azimuth, and analysis point are smoothed by performing a running window time 

integration over the partial sums from -K to -^K: 

-kkV J > 
k=-KV j=l 

+lc J " 

2!SK0 + kAt,p, q, Xj,yj)f 

k=.K M 

where K is the half width of the t!me window in samples. 

10. The method of Claim 1. wherein said traces within said analysis cells ar 
15 characterized by a maximum dip and a maximum temporal frequency component; and 

wherein step (c) includes the steps of: 

obtaining an estimate of the maximum true dip and the maximum temporal 
frequency component of said traces in said analysis cell; 



RECTIFIED SHEET (RULE 91) 
ISA/EP 



BNSDOCID: <WO_9713166A1J_> 



wo 97/13166 



• 34- 



PCTAJS96/15705 



using said maximum true dip, said maximum temporal frequency and said 
pre-determined lateral dimensions of said analysis cell to calculate apparent dip 
increments in two generally perpendicular directions relative to said dip azimuth 
measurement axis. 

5 11. The method of Claim 1 . where in performing step (c) said measure is at 

least a function of: 



where J Is the number of traces in said analysis cell, where Uj(T, p. q) is a representation 
of the seismic trace in said analysis cell, where t is the time, p is the apparent dip in the x 
10 direction, and q is the apparent dip in the y direction; wherein p and q are measured in 
ms/m and the x and y directions are mutually perpendicular. 

12. The method of Claim 1 1 . where in performing step (c) said measure is also 
a function of. 

ZMt-(px,+^,)]}^ 

"•5 13. The method of Claim 12. where in perfbmiing step (c) said measure is a 

function of: 
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14. A method of locating subterranean f atures, faults, and contours, 
comprising the steps of: 

(a) accessing 3D seismic data covering a pre-determined volume of 

the earth; 

5 (b) dividing said volume into an array of relatively small three- 

dimensional cells wherein each of said cells is characterized by at least five lat rally 
separated and generally vertical seismic traces located therein; 

(c) determining in each of said cells the semblance/similarity of said 
traces relative to two pre-determined directions; and 

10 (d) recording said semblance/similarity of said cells in a form for 

display as a two-dimensional map of subterranean features. 

15. The method of Claim 14. where in performing step (c) said pre-determined 
directions are mutually perpendicular; and said semblance/similarity of said traces within 
each cell is a function of at least time, the numt)er of seismic traces within said analysis 

15 cell, and the apparent dip and apparent dip azimuth of said traces within said analysis 
ceil. 

16. The method of Claim 15, where said semblance/similarity of said traces 
within each cell is determined by computing a plurality of measurements of the 
semblance/similarity of said traces within each cell and selecting the largest of said 

20 measurements of said semblance/similarity of each cell; and wherein step (c) further 
includes the step of defining the apparent dip and apparent dip azimuth corresponding to 
said largest of said measurements to be an estimate of the true dip and an estimate of 
the true dip azimuth of the seismic traces within said analysis cell. 
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17. The method of Claim 16, wherein each of said plurality of measur ments 
of said semblanc /similarity is at least a function of the en rgy of said trac s; and wher in 
said energy of said traces is a function of time, the number of seismic traces within said 
analysis cell, and the apparent dip and apparent dip azimuth of said traces within said 

5 analysis cell. 

18. The method of Claim 16. wherein said map is a color map that is 
characterized by hue, saturation and lightness; 

wherein one of said estimates of true dip azimuth, said estimates of true dip, and 
said largest calculated measures of semblance is mapped onto one of a lightness scale. 
10 hue scale, and a saturation scale; 

wherein another of said estimates of true dip azimuth, said estimates of true dip. 
and said largest calculated measures of semblance is mapped onto another of said 
lightness scale, said hue scale, and said saturation scale; and 

wherein the remaining one of said estimates of true dip azimuth, said estimates of 
15 true dip, and said largest calculated measures of semblance is mapped onto th 
remaining one of said lightness scale, said hue scale, and said saturation scale. 

19. The method of Claim 18, wherein step (d) comprises the steps f: 
mapping said estimates of true dip azimuth onto said hue scale, mapping said estimat s 
of true dip onto said saturation scale, and mapping said largest calculated measures of 

20 semblance onto a lightness scale. 

20. In seismic exploration wherein 3D seismic data comprising reflected 
seismic energy is recorded as a function of time and wherein a computer is used that is 
programmed to process such seismic traces and to produce an image therefrom that is 
representative of subterranean features, an article of manufacture comprising: 
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a medium that is readable by a computer and that carries instructions for said 
computer to perform a process comprising the steps of: 

(a) accessing 3D seismic data over a predetermined volume of the 
earth, said data comprising seismic traces that are characterized by time, position and 

5 amplitude; and 

(b) ascertaining the similarity of nearby regions of said 3D seismic data 
of said volume by: 

(1) dividing at least a portion of said data into an array of 
relatively small, adjacent, three-dimensional analysis cells, wherein each of said analysis 

10 cells contains portions of at least five seismic traces; and 

(2) computing a seismic attribute for each cell that is a function 
of the largest of a plurality of measurements of semblance and the corresponding 
apparent dtp and the corresponding apparent dip azimuth. 

21. The article of manufacture of claim 20, wherein said medium carri s 
IS instructions for the computer to perform step (2) by making measurements of semblance 
that are a function of: 

J 

where x and y are distances measured from the center of the analysis cell along mutually 
perpendicular x and y axes, where J traces is the number of seismic traces, where u^(x, 
20 p. q) represents a seismic trace, where t is the time, p is the apparent dip in the x 
direction, and q is the apparent dip in the y direction; and wherein p and q are measur d 
in ms/meter. 
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22. The article of manufacture of claim 
instructions for the computer to perform step (2) 
semblance that are also a function of : 

J 

5 23. The article of manufacture of claim 21, wherein said medium carries 

instructions for said computer to perform step (1) by forming analysis cells having an 
elliptical cross-section. 

24. The article of manufacture of claim 23. wherein said predetermined volum 
is characterized by a fracture having an ascertainable direction: and wherein said medium 

10 carries instructions for said computer to form analysis ceils that are generally elliptical in 
shape and that have major axes aligned in the direction of said fracture. 

25. In seismic exploration wherein reflected seismic energy is recorded as a 
function of time to produce a series of seismic traces, a method comprising the steps of: 

(a) accessing a data set of seismic traces distributed over a three* 
15 dimensional volume of the earth, said volume of the earth having subterranean features 

characterized by dip and dip azimuth: 

(b) calculating a plurality of measures of the semblance of said traces 
within a relatively small three dimensional analysis cell that is located within said volume 
and at one part of a predetermined time layer, wherein each measure of semblance is at 

20 least a function of time, the number of seismic traces within said analysis cell, and the 
apparent dip and apparent dip azimuth of said traces within said analysis cell; 

(c) computing a seismic attribute for said analysis cell that is at least a 
function of the largest of said plurality of calculated measures of semblance and the 
corresponding apparent dip and the corresponding apparent dip azimuth, wherein said 

25 corresponding apparent dip and said corresponding apparent dip azimuth are defined to 



PCTAJS96/15705 

21. wherein said m dium carries 
by making measurements of the 
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be stimates of th true dip and an stimate of the true dip azimuth of the seismic traced 
within said analysis cell; 

(d) repeating steps (b) and (c) along other parts of said time layer, and 

(e) forming a map of said seismic attributes over said time layer. 
5 26. The method of Claim 25. wherein step (a) comprises the steps of: 

(1) accessing 3D seismic data over a predetermined volume of the 
earth, said 3D seismic data comprising at least eleven seismic traces that are 
characterized by time, position and amplitude: and 

(2) dividing a portion of said volume into at least one time lay r 
10 comprising an an^y of relatively small, three-dimensional cubes that contain at least five 

seismic traces; and wherein said cubes are used as the cells to perfomi step (b). 

27. The method of Claim 26. where in perfomiing step (b) each measure of 
semblance is a function of: 

j=' 

—J 

15 where each analysis cell contains portions of at least J seismic traces, where J is at least 
5. where x and y are distances measured from the center of the analysis cell along 
mutually perpendicular x and y axes, where p and q are the apparent dips in the x and y 
directions, where u^t. p. q. x. y) represents a seismic trace within said analysis cell, and 
where the true dip d and dip azimuth tp are related to p and q by p = d sin(q>) and q = d 

20 cos(<{>). 

28. The method of Claim 27. wherein each measure of semblance for each 
dip. dip azimuth, and analysis point are smoothed by fonning a njnning window time 
integration over partial sums of a time window within said horizontal time layer. 
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29. A method of seismic exploration, comprising the steps of: 

(a) r ading a 3D seismic data set comprising seismic signal traces that 
are distributed over a volume of the earth; 

(b) selecting at least one horizon slice from said volume and forming 
5 therein cells that are arranged into laterally extending rows and columns, each of said 

cells having at least five seismic traces extending generally therethrough; 

(c) computing for each of said cells: 

(1) a plurality of semblance measurements of said traces, 
wherein each measurement is at least a function of time, the number of seismic trac s 

10 within said analysis cell, and the apparent dip and apparent dip azimuth of said traces; 

(2) the largest of said plurality of measurements of semblance; 

and 

(3) an estimate of the true dip and an estimate of the true dip 
azimuth of the seismic traces within said analysis cell from the apparent dip and apparent 

15 dip azimuth corresponding to said largest measurement of semblance; and 

(d) displaying, over said at least one horizon slice, of representations 
of said largest measurements of semblance and said estimated true dips and said 
estimated true dip azimuths of each of said cells. 

30. The method of Claim 29. wherein step (b) is performed by selecting a 
20 horizon slice that is characterized by a common time; and wherein step (d) is performed 

by displaying across said time slice representations of said largest measurements of 
semblance and said estimated true dips and said estimated true dip azimuths of said 
cells. 

31. The method of Claim 29. wherein step (d) is performed by forming a color 
25 map that is characterized by hue. saturation and lightness, wherein for e?ch of said cells: 
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one of said stimat s of true dip azimuth, said stimates of true dip. and said 
larg st calculated measurements of semblanc is mapped onto one of a lightness scale, 
hue scale, and a saturation scale; 

wherein another of said estimates of true dip azimuth, said estimates of true dip. 
5 and said largest calculated measurements of semblance is mapped onto another of said 
lightness scale, said hue scale, and said saturation scale; and 

wherein the remaining one of said estimates of true dip azimuth, said estimat s of 
true dip. and said largest calculated measurements of semblance is mapped onto the 
remaining one of said lightness scale, said hue scale, and said saturation scale. 

10 32. In the exploration for gas and oil wherein over a volume of the earth 

seismic traces are recorded, a method comprising the steps of: 

(a) grouping at least parts of at least five relatively close seismic trac s 
into a plurality of relatively small three-dimensional analysis cells; 

(b) performing in each of said cells a plurality of measurements of the 
15 semblance of said parts of said traces as a function of at least time, the number trac s 

therein, the apparent dip of said traces, and the apparent dtp azimuth; 

(c) identifying in each of said cells the largest of said plurality of 
measurements of semblance, the corresponding apparent dip, and the corresponding dip 
azimuth: and 

20 (d) converting said largest measurements of semblance, said 

corresponding dip and said corresponding dip azimuth of said cells into color attributes of 
hue. saturation and lightness, wherein for each cell: 

one of said dip azimuth, said dip, and said largest measurements of 
semblance is mapped onto one of a lightness scale, hue scale, and a saturation scale; 
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another of said dip azimuth, said dip. and said largest m asurements of 
semblance is mapped onto another of said lightness scale, said hue scale, and said 
saturation scale; and 

the remaining one of said dip azimuth, said dip, and said largest 
5 measurements of semblance is mapped onto the remaining one of said lightness scale, 
said hue scale, and said saturation scale. 

33. A device adapted for use by a workstation wherein 3D seismic data is r ad 
into memory and processed into a color display of subterranean features, comprising: » 

computer readable means carrying instructions for a process comprising the steps 

10 of: 

(1) digitally locating said 3D seismic data in an array of relatively small 
three-dimensional cells, wherein each of said cells contains representations of a part of at 
least five seismic traces; 

(2) calculating for each of said cells an estimate of the semblance, an 
15 estimate of the true dip. and an estimate of the true dip azimuth of said parts; and 

(3) converting said estimates of semblance, said estimates of true dip. 
and said estimates of true dip azimuth into an array of digital values corresponding to the 
color attributes of hue. saturation, and lightness, 

34. The device of Claim 33. wherein one of said estimates of true dip azimuth. 
20 said estimates of true dip. and said estimates of semblance is mapped onto one of a 

lightness scale, a hue scale, and a saturation scale for each of said cells: 

wherein another of said estimates of true dip azimuth, said estimates of true dip. 
and said estimates of semblance is mapped onto another of said lightness scale, said hue 
scale, and said saturation scale for each of said cells; and 
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wherein the remaining one of said estimates of true dip azimuth, said estimates of 
true dip. and said estimates of semblance is mapped onto the remaining one of said 
lightness scale, said hue scale, and said saturation scale for each of said cells. 

35. The device of Claim 33. wherein said computer readable means carri s 
5 instructions to perform step (2) by: 

(1) calculating a plurality of semblance measurements relative to at least 
two directions, and selecting the largest of said measurements; 

(ii) selecting the apparent dip corresponding to said largest measur m nt 
of semblance from step (i); and 

10 (iii) selecting the apparent dip azimuth corresponding to said largest 

measurement of semblance from step (i). 

36. The device of Claim 33. wherein said computer-readable means is sel ct d 
from the group consisting of a magnetic tape, a magnetic disk, an optical disk and a 
CD-ROM. 

15 37. A method of prospecting for hydrocarbon deposits, comprising the steps of: 

(a) obtaining a color seismic attribute display of 3D seismic data for a pre- 
determined three-dimensional volume of the earth, said display being generated by using 
data obtained by a computer and at least one program for said computer that instructs said 
computer to perform the following steps: 

20 (1 ) convert said volume into an array of relatively small three-dimensional 

cells, wherein each of said cells has a portion of at least five seismic traces located therein; 

(2) make a plurality of semblance measurements within each of said 
cells, wherein each measurement is at least a function of time, the number of seismic 
traces within said cell, the apparent dip of said traces and apparent dip azimuth of said 

25 traces; 
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(3) select the largest of said plurality of measurements of semblance in 

each cell; 

(4) use as an estimate of the true dip and an estimate of the true dip 
azimuth in each cell the apparent dip and apparent dip azimuth that correspond to said 

5 largest measurement of semblance in said cell; 

(5) map said estimates of true dip azimuth onto a hue scale; 

(6) map said estimates of true dip onto a saturation scale; and 

(7) map said largest calculated measures of semblance onto a 
lightness scale; and 

10 (b) using said color display to identify subsurface structural and sedimentologtcal 

features commonly associated with the entrapment and storage of hydrocarbons. 

38. The method of Claim 37, further including the step of using said map to 
identify drilling hazards. 

39. The method of Claim 38. further including the step of drilling at a location 
15 identified in step (b). 

40. The method of Clairri 37, wherein step (a)(2) comprises the step of 
computing: 

J 

X • (px j + qy j)] 

where each cell is characterized by two perpendicular dimensions, where x and y are 
20 distances measured from the center of the cell along mutually perpendicular x and y 
axes, where J is the number of seismic traces, where Uj(T, p. q) represents a seismic 
trace, where x is the time, p is the apparent dip in the x direction, and where q is the 
apparent dip in the y direction. 
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41. The method of Claim 40, wherein step (a)(2) comprises the step of 
computing: 

42. A map for the exploration of oil and gas, comprising: 

5 (a) a generally flat medium for recording visually perceptible color images 

thereon, said color images being characterized by hue, saturation and lightness; and 

(b) a plurality of images on said medium that are representative of semblance, 
true dip and true dip azimuth of 3D seismic data from a pre-determined volume of the earth, 
wherein true dip azimuth is mapped onto a hue scale, wherein true dip is mapped onto a 
10 saturation scale, and semblance is mapped onto a lightness scale. 

43. The map of Claim 42. wherein said volume Is representative of an areal 
array of three dimensional cells; 

wherein each of said cells contains a representation of parts of at least fiv 
seismic traces; 

15 wherein each of said images correspond to one of said cells; 

wherein the semblance of the traces within each cell is determined by making a 
plurality of measurements of the semblance of said traces within each cell and selecting 
the largest of said measurements; 

wherein each measurement of semblance is at least a function of time, the 
20 number of seismic traces within said cell, and the apparent dip and apparent dip azimuth 
of said traces: 

wherein the true dip of each cell is represented by the apparent dip that 
corresponds to said largest of said measurements of semblance; and 
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wherein th tru dip azimuth is represented by the apparent dip azimuth that 
corresponds to said iarg st of said measurements of semblance. 

44. The map of Claim 43. wherein semblance is measured as a function of: 

J 

uf-f - (px j + qy j)] 

5 where x and y are distances measured from the center of each cell along mutually 
perpendicular x and y axes, where J is the number of seismic traces located ther in. 
MUhere Uj(x. p. q) represents a seismic trace, where t is the time, p is the apparent dip in 
the X direction, and q is the apparent dip in the y direction. 

45. The map of Claim 44. wherein semblance is measured as a function of: 

{Z " - (px j + qy j)] y 

j-i 



10 



46. The map of Claim 42. wherein said medium is the face of a cathode ray tube. 

47 In a computer workstation wherein 3-D seismic data obtained over a pre- 
determined three-dimensional volume of the earth is read into n^emory. wherein a computer 

15 divides such volume into an array of three-dimensional analysis cells, wherein each cell has 
at least a portion of five laterally separated seismic traces located therein, and wherein the 
computer is used to transfomi such data into a display of seismic attributes, the computer 
CHARACTERIZED BY performing a process comprising the steps of: 

(1) calculating in each of the cells a semblance value for said seismic traces. 

20 wherein said semblance value is at least a function of time, the number of seismic trac s 
within said cell, the apparent dip of said traces, and the apparent dip azimuth of said traces: 
and 
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(2) displaying said semblance value of each cell that lies between two planes 
within the 3-D volume to identify subsurface features commonly associated with the 
entrapment and storage of hydrocarbons. 

48. The computer workstation of Claim 47, wherein the computer performs step 
5 (1) by: making a plurality of semblance measurements within each of said cells; and 

selecting the largest of said plurality of measurements as said semblance value of said cell. 

49. The computer workstation of Claim 48. wherein after performing step (1) the 
computer performs the step of: using the apparent dip and the apparent dip azimuth that 
correspond to said largest measurement of semblance in said cell as an estimate of tru dip 

10 and as an estimate of true dip azimuth of said ceil. 

50. The computer workstation of claim 49. wherein the display of step (2) is 
characterized by color components of hue. saturation and lightness; and wherein step (2) 
comprises the steps of mapping said estimate of true dip azimuth for each cell onto a hue 
scale; mapping said estimate of true dip for each cell onto a saturation scale; and 

15 mapping said largest calculated measures of semblance onto a lightness scale. 
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